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Abstract 

Two-pion exchange parity-violating nucleon-nucleon interactions from recent effec- 
tive field theories and earlier fully covariant approaches are investigated. The po- 
tentials are compared with the idea to obtain better insight on the role of low-energy 
constants appearing in the effective field theory approach and the convergence of 
this one in terms of a perturbative series. The results are illustrated by considering 
the longitudinal asymmetry of polarized protons scattering off protons, p+p p+p, 
and the asymmetry of the photon emission in radiative capture of polarized neutrons 
by protons, n + p ^ d + ^. 



1 Introduction 



Effective field theory (EFT) underlies most recent developments in the domain of the 
nucleon-nucleon [NN) strong interaction [H [2l El H] . The approach is mainly motivated by 
the fact that a large part of the short-range interaction is essentially unknown. Its detailed 
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description may not be really relevant at low energy and a schematic one, represented 
by contact interactions with low-energy constants (LEC's), could be sufficient. Moreover, 
such an approach could account for important properties in relation to QCD dynamics, i.e. 
chiral symmetry. Implementing these properties can be done with the chiral perturbation 
theory [S] . One is thus led to distinguish contributions at different orders. Beyond the one- 
pion exchange (OPE), which appears at leading order (LO), two-pion exchange (TPE), 
which is relevant at next-to-next-to- leading order (NNLO), has been considered^ Higher 
order terms are also considered, contributing to a successful description of the strong 
interaction. 

Naturally, the EFT approach has been applied to the weak, parity- violating (PV), 
NN interaction [6J, superposing on earlier phenomenological works in the 70's [3 [8], [9] 
with a systematic perturbation scheme in terms of an expansion parameter characterizing 
the theory. Thus, a one-pion-exchange contribution appears at LO while the two-pion- 
exchange contribution is part of those at NNLO. It was claimed that effects from the 
two-pion-exchange contribution could be potentially large. Estimates have been made 
for various observables [TOl [11] and, while they do not contradict the above expectation, 
they have evidenced a large range of uncertainty ^U\. This points out to the role of 
a contact term present in the operators, which has to be completed in any case by a 
LEG contribution and cannot be therefore considered as physically relevant by itself. 
The PV TPE contribution was considered in the 70's in several works starting from a 
covariant formalism and based on Feynman diagrams [12] or dispersion relations [131 [El 
[15] . Originally, these works were motivated by the expectation that the TPE could play 
a role in the PV case as important as in the strong interaction one, what was actually 
disproved by the studies. A similar motivation has recently been addressed within the 
EFT approach with some attention to the contribution involving the A excitation [16j 
(see also Refs. [TTl [18] with this last respect). Interestingly, the TPE contribution in 
various processes turned out to be well determined but rather small [13l [15l [19] and, in 
particular, unessential in comparison with other uncertainties (PV couplings constants, 
nuclear- structure description, etc). This last feature largely explains its omission in later 
works. In view of different conclusions, we believe that a comparison of both the recent 
and earlier works is useful. Some preliminary results were presented in Ref. |2QJ . 

The above studies could be relevant because, contrary to the strong interaction, it is 
not possible to determine at present the LEC's due to the lack of sufficient and accurate 
enough experimental data. On the one hand, they could tell us about the role of contact 
interactions in making the two approaches as close as possible. The contributions of these 
contact interactions can be ascribed to the LEG or the "finite" range part, depending 
on the subtraction scheme. On the other hand, they could provide information on the 
convergence of the perturbative expansion of the potential in the EFT approach, which 
is limited to NNLO so far. In the field of the strong NN interaction [21] or weak semi- 
leptonic interactions [22], there are hints for non-negligible corrections. 

The plan of this paper is as follows. In the second section, results relative to the TPE 
contribution from a covariant approach are reminded (isovector part). All components 

^ There are different conventions to denote orders. We use the one in agreement with what has been 
used in the parity-violating case. 



2 



of this interaction and their convergence properties are in particular discussed. The 
third section is devoted to results in the EFT approach, completed by those obtained 
from the contribution of time-ordered diagrams as a check. In the fourth section, we 
examine similarities and differences between results of these approaches and those from 
an expansion of the covariant one at the lowest non-zero order in the inverse of the nucleon 
mass, 1/M. A numerical comparison concerning a few aspects of potentials so obtained is 
given in the fifth section. Estimates of the effects in two selected processes, proton-proton 
scattering and radiative neutron-proton capture at thermal energy, are finally given in 
the sixth section. These two processes allow one to illustrate the two types of PV effects 
that are expected from the TPE interaction at low energy. The seventh section contains 
the conclusion. This is completed by appendices concerning the removing of the iterated 
one-pion exchange and the EFT approach. 



2 Two-pion exchange from a covariant formalism 



X. 





Figure 1: Two-pion exchange in the covariant approach: These diagrams represent the 
contributions of the crossed box, the non-crossed box and the iterated OPE that has to be 
subtracted from the previous one. The continuous fine represents a baryon and the dashed 
one a pion. The contributions with nucleon and nucleon resonances in the intermediate 
state have been considered in the literature. Only the first one is retained here but the role 
of the other ones will also be mentioned. The PV meson-nucleon vertex is marked with 
a cross, x. The last diagram on the right represents the iterated one-pion exchange, the 
back slash indicating that the corresponding nucleon is on-mass shell. Further diagrams 
with a different order of the parity- violating and parity-conserving meson-nucleon vertices 
on the same nucleon line, not shown here, are also considered. 

We here consider the isovector PV TPE contribution to the NN interaction obtained 
from a fully covariant formalism. It is induced by an elementary PV pion-nucleon coupling, 
most often denoted by h^^. As this coupling also determines the strength of the PV OPE 
interaction, we give the corresponding expression in terms of both the momentum transfer, 
q, and the relative momenta, p and p': 

Due to possible ambiguity, we further specify the notations for momenta: 

1 1 

Q = {Pi- P'l) = -{P2- p'2) = {P- P') , (2) 
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where the primed and non-primed momenta refer to those of particles appearing respec- 
tively in the bra and ket states. In configuration space, the interaction thus recovers its 
standard form: 



1 Q-kNN 



hi 



47rr 



(ri X ((Ti + a2)-(ri-r2) ^ , (3) 



2V2M ^' ' ' ^' 47rr3 

with r =\f\— f2\. First studies of the isovector PV TPE contribution were made in the 
early 70's [HI |13l [HI |15]. For a part, they were motivated by the underestimation of 
some observed PV effects using the standard PV A^A^ interaction available at that time. 
Similarly to the strong-interaction case, where the TPE contribution is quantitatively 
more important than the OPE one, it was believed that the TPE contribution could also 
play an essential role in the PV case. 

Various studies roughly agree with each other, after correcting mistakes in some cases 
[121 [SI . Differences involve in particular the formalism (calculations using Feynman dia- 
grams or dispersion relations), non-relativistic approximations in external nucleon lines, 
or the removing of the OPE-iterated contribution in the box diagram. The choice of the 
Green's function in the last ingredient was rather unimportant due to the introduction 
of cutoffs in applications. It could however be important in unrestricted calculations. 
The point is of relevance with respect to a comment made in Ref. [B] about the absence 
of convergence in earlier calculations. It will be discussed in more detail below when 
expressions for the TPE contribution are given. 

The crossed and non-crossed box diagrams that enter the isovector TPE contribution 
of interest here are represented in Fig. [H where the intermediate hadron on the upper 
line can be a nucleon as well as a resonance (A(1231), A^*(1440), A^*(1518)). For our 
purpose and also for simplicity, we only retain the nucleon. Altogether, the corresponding 
contribution to the isovector interaction involves six different terms. Following for a part 
the notations of Ref. [12], it can be written quite generally in momentum space as: 



y{v\ P) = Vu + Vu + ^56 + ^75 + ^66 + ^85 

= i (ri + T2Y {ai X 02) ■ {p - p) vu{q, ■ ■ ■) 
+ (ri + T2y (ai - 0^2) -{p +P) V34{q, ■ ■ ■) 
+i (ri X T2y (ai + ^2) -{p' -P) v^eiq, ■ ■ ■) 
+ in - T2y (ai + (T2) -{p +p) V75{q, ■ ■ ■) 
+ (rixr2)^ (ffi-g CT2-(^'+p) xg + (ai ^ (T2)) Veeiq,--) 

-i{Ti~T2Y [ai-{p'+p) a2-{p'+p)xq+ {ai ^ a2)) vs^iq,- ■■) ■ (4) 
Functions Vij{q, ■ ■ ■) assume a dispersion-relation form: 

■' TT J4m2 (t' + g2) 

and dots represent possible extra dependence on p' and p (kinetic energy in particular). 
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The configuration-space PV TPE potential is obtained from tlie standard relation: 



v{r) = I (6) 

Its expression tlius reads: 

V{r, p, p) = i (ri + (<?i X 0^2) ■ [p, Viiir, ■ ■ ■)] 
+ + ■r2y (ai - 0^2) -{p, v^iir, ■ ■ ■)} 
+i (ri X (ai + 0^2) ■ [p, t^56(r, ■ ■ ■)] 

+ - -r2y (cti + (T2)-{P, ^^75(^, ■ ■ ■)} 

+2i {TiyiT2y {Br[p, ^2-1 -—veeir, ■■■)] + {^1 ^ 0^2)) 
\ r ar ' 

-2 (ri-ra)" ^2 • r ^ -^t;85 (r, ■■■)} + (^^i ^ '?2)) , (7) 

^ r ar ' 

where: 

while p' and p, which dots account for, have now an operator character and should be 
placed respectively on the left and the right, in accordance with our conventions. 

The terms V44, V34, V56 and V75 appear at the first order in a p/M expansion of the 
Lorentz invariants appearing in the expression of the interaction. The two other ones, 
Vge and Vgs, appear at the third order in p/M. They can therefore be considered as 
relativistic corrections. Moreover, they only contribute when going beyond the transitions 
between lowest partial wave states, i.e. S to P, which generally dominate at low energy 
(where most PV data are available). For these two reasons, the corresponding terms were 
discarded in the past, which we also do here. We however stress that these higher-order 
terms are necessary to get a full mapping of the NN interaction, especially to discriminate 
transitions involving higher partial waves such as ^Pi — ^Di and ^P2~^-D2, ^-D2~^-^2 and 
^Ds — ^Fa, etc. It is also noticed that the above non-relativistic expansion only involves 
the nucleon external lines, as done in other approaches. None is made for internal nucleon 
lines where big effects could arise. As our results presented in this paper only retain part 
of the full relativistic structure, they will be denoted "covariant" to avoid overstating this 
property. 

We now consider the expressions of the V44, V34, V56 and 1^75 terms. It is first noticed 
that V44 and V34 only receive contribution from the crossed diagram while V56 and V75 
also get some from the non-crossed one. In this case, one has therefore to worry about the 
removing of the iterated OPE and, especially, about the choice of the Green's function 
which enters the calculation. In Ref. [13], the Green's function, {2Eq — 2E)~^, was used, 
taking into account that the Schrodinger equation is linear in the energy of the system, 
Eq, and assuming moreover that the kinetic energy of particles retains its relativistic 
form {E = y/ M'^ + p'^ in the cm.). In later works [HI [TB], a Green's function more in 
agreement with a non-relativistic Schrodinger equation, E/{pl — p^), was instead used. 
The difference, a factor 2E/[Eo + E), had minor numerical effects in the past calculations 
where cutoffs were introduced in the dispersion integrals, Eq. ([5]). The difference however 
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matters in unrestricted calculations. Dispersion-relation integrals diverge in the first case 
while they converge in the second one. With this last respect, it is noticed that the 
corresponding Green's function can be written in the following form, E/{Eq — E"^), which 
rather evokes an equation with a quadratic-mass operator. Such an equation is known to 
provide solutions with a behavior in the relativistic domain better than the linear one [23j . 
This is the choice made here. Taking into account that the kinetic-energy dependence 
is small, and that we are interested in low-energy processes, the momentum in functions 
g{q, ■ ■ ■) can be set to 0. In a similar way as Ref. |13] with details given in Appendix lA.ll 
the closed expressions for these functions are then obtained. Omitting dots that are no 
longer justified, they read: 




H X 



2 



2 + 4M2g2 



X 



~ 1 H X X / 77^2 \ 

556(0 = -K— ^ ^^^^^ - ^'^H^J 



'kH'-{ml + k^y 2E(E+M)\M^ M{E+M) 



r> / ml ^ ~ 2G 



-K 



■lkH'-{ml + k^Y 



2 / x2 2Ex _ e _ k^ \ 
^ E{E + M) \MH'^ Mf ~ M{E+M)~ M{E + M)t'j ' ^' 



where: 



K 



3271^2 




^ ^ 2 / ^' + 4M2g2 / ^x^ + AM^q^^ + q^Vf \ 
V t' \^x^+AM^ql-q^y/V) ' 

kl = x±q^Vt', E = VWT¥. (10) 



Looking at the asymptotic behavior of the functions gij{t') for large t', it is noticed that 
the dominant contributions of individual terms, oc t'^^'^, as well as the constant ones, cancel 
(see Appendix IA.2I for some detail). One is thus left with the following contributions: 

^44(0.^00 = K^(ln(^)-l), 
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1 /5 / t' X 15 3 , , A 
= -*M7Fl8"<M^)+T-2'"P)j' 

= /^_L^g,„(^)-l + to(2))^ (11) 

In all cases, the functions gij{t') behave asymptotically as up to log terms, ensuring 

the convergence of the integrals given in Eq. ([5]). This result is important as it allows one 
to consider the dispersion approach as a benchmark, thus providing information about 
contributions that are ascribed to LEC's as well as possible higher-order corrections in 
other approaches. 



3 Two-pion exchange from the EFT approach 



3.1 EFT approach 




Figure 2: One and two-pion exchange in the EFT approach: The contributions from the 
OPE (a), from the triangle TPE (b), from the crossed TPE (c), from the non-crossed 
TPE (d) and from the contact term (e). See the caption of Fig. [T]for further comments. 

The PV TPE at NNLO has been calculated in the EFT approach by Zhu et al. [6]. It 
contains two components which, in their notations, involve functions C*|'^(q') and C|'^(g) 
and correspond here to the components V44 and V^q of the more complete interaction given 
by Eq. (jl]). The contributions being accounted for correspond to the diagrams shown in 
Fig. [21 Their expressions for the finite-range part and the associated contact term have 
been obtained in the maximal-subtraction (MX) scheme. Factorizing out the spin-isospin 
dependence and taking into account corrections made since then [TU] , they read: d 

= -4v/2vr^(,iL(g)), 

X 

vT{<l) = -V2 7.^[gA m~g\ {3L{q)-H{q))j , (12) 

where the scale is roughly given by = 47r/^ ~ An qaM / Qt^jsin — 1 GeV. The L[q) 

^We have followed our own conventions, as the conventions used to present the final results in Ref. [5] 
differ from that ones defined by the same authors at an earlier stage. 
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and H{q) functions are defined as: 



^q'^ + Aml ^ / yJq'^+Aml + q\ 
2? \^q^+4ml - q) 

(13) 

The detail of the contributions corresponding to the diagrams (b), (c) and (d) in Fig. [2]can 
be found in Appendix [Bl The above terms entering the interaction should be completed 
by contact terms: 

v'iI = C,,. (14) 

The contributions from two-pion exchange to these contact terms are also given in Ap- 
pendix [HI where it is seen that they require some renormalization. The sum of the EFT 
two-pion-exchange and contact contributions should be well determined, but how it is 
split between the two terms is not. In the minimal-subtraction (MS) schemejl the part of 
the contact term which is proportional to / 711^^) (see Appendix [B|) is shifted to the 

EFT two-pion-exchange part. In such a case, the term "1" cancels the function L{q) at 
g = and the log term, for yU > tjIt^, changes the overall sign of the potential in the low-g 
range. 




3.2 Relation to the time-ordered-diagram approach 

When integrating the EFT expressions of the TPE interaction over the time component 
of the integration variable entering some loop, it is expected that one should recover 
expressions obtained from considering time-ordered (TO) diagrams in the non-relativistic 
limit ({7 — > where v is the nucleon velocity v ~ p/M). As this check was quite useful in 
determining the correct expressions given in the previous subsection, we give below the 
raw expressions obtained from the contribution of these time-ordered diagrams, which are 
shown in Fig. [31 Starting from the elementary nNN interaction, q^nn^ ■ k/{2M y/uJk) 



(with Uk = \Jml. + A;2), one gets: 
y-nNN "-TT 



4^2 M3 

dk k^-(k-qf 



+ 7 — : ^ + 



4°(^) 



{2tiY OJiOJj \Ujf [uJi + UJj) {uJi + UJj) Uj'j Ui (uJi + UJj) Uj 

qInnK fl f dk 1 



4^2 M3 {27r)^UiUJj{uJi+ujj) 
1 r dk A;2-4 / 1 



2 J {2tt)^ ujiUJj \uj]; {(jJi+ujj) {uji+ujj) u'j uJi{ijJi+ujj) uj 



, (15) 



■^What we denoted by MN the scheme in a previous work [10] actuaUy corresponds to the MS scheme, 
hence the change of name adopted here. 
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Figure 3: Two-pion exchange in the time-ordered diagram approach: The contributions 
from the crossed diagrams (the first hne), from the Z-type contributions (the second hne, 
identified as triangle diagrams in a different approach), from the non-crossed diagrams 
(the diagrams on the third hne). The last two diagrams differ in that the first of them 
involves a meson propagator with off-energy shell contributions, which are omitted in the 
other one (what is reminded by a back slash on the meson lines). This last contribution, 
which arises from the iterated OPE, has to be subtracted from the previous one. The 
discrepancy involves a factor, Eq — k'^/M, which cancels a similar factor entering the 
denominator of the Green function {k can be identified as the loop momentum). See the 
caption of Fig. [T] for further comments. 



where ooi = ym'^ + {k + Uj = ym'^ + {k — q/2y. The three integrals involve the 

contributions successively from the crossed diagrams (the first line of Fig. [3]), Z-type ones 
(the second line of Fig. [3]) and both crossed and non-crossed diagrams (the first and third 
lines of Fig. [3]). The Z-type diagrams are calculated assuming a pseudo-scalar coupling, 
consistently with the dispersion-relation approach used independently, and retaining the 
lowest non-zero term in a 1/M expansion. As is known, the contribution alone violates 
chiral symmetry (see, for instance, Ref. [21]). The expected symmetry is restored by 
further contributions, which can be calculated in the same formalism (see some detail in 
Sec. 14.21) . Contributions of all diagrams in Fig. [3] diverge but they contain a well-defined 
part that can be analytically calculated. Interestingly, the expressions so obtained can 
be cast into the form of dispersion integrals. This property stems from considering a 
complete set of topologically-equivalent time-ordered diagrams. This writing is interesting 
as it greatly facilitates the comparison with the expressions obtained from a covariant 
approach, which evidences the same form. It is thus found that the different integrals in 
Eqs. f|T5l) read: 

k^-{k-qf (I 1 1 \ 

dk ^ ^ + — . + ' 



An(l~L(q))+ dk^ = 7r dt' 



dk 
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27r(l - L(q)) + - dk^ = TT dt' 



I 



dk ^ ^r- + + 



UJi UJj \ Ujf {uJi + UJj) {uJi + UJj )uj'j UJi{uJi + UJj ) Uj ^ 

27r [3 (l -LW) +//(,))+ II <i*^ 



J4n 



df ^(t'-An^D+Anj ^^^^ 



^mi ^ Jt'-Amlit' + q^) 



Due to the divergent character of these integrals, the above equahties hold up to some 
constant. This does not however affect the g^-dependent part. One can thus remove an 
infinite contribution so that the interaction takes a definite value at some q^. A particular 
interesting choice is to subtract a part so that the remaining one, which contains the most 
physically relevant part, cancels at = 0. By looking at this quantity, one can usefully 
compare different approaches. To some extent the slope with respect to at = 
provides information on the sign and the strength of the interaction at finite distances. 
The L{q) function in the above equalities thus points to a configuration-space interaction 
with an opposite sign at finite distances. The divergent part is not without interest 
however. It tells us in which direction the (short-range) subtracted interaction is likely to 
contribute. By integrating out the contribution t' > h? at the r.h.s. of Eqs. (fT6|) in the 
small q limit, one successively gets the approximate factors 47r ln(A/2m^), 27r ln(A/2m^), 
and 67rln(A/2m7r). This suggests two observations. On the one hand, for large enough A, 
the r.h.s. has a sign opposite to that one given by the L(g) term at the l.h.s., confirming 
the observation in configuration space. On the other hand, the above factor allows one 
to make some relation with the EFT interaction calculated in the minimal-subtraction 
scheme, which involves similar log terms (with A replaced by yu). 

Comparing the above expressions, Eqs. f|T6l) . with the previous EFT ones, Eqs. f lT2l) . it 
is found that the q^ dependences are very similar. Assuming the choice A^ = 47r qa M/ Qt^nn, 
an identity is actually found for the potential, f 44(g), as well as the second term for the 
other potential, v^Q^q). For the first term in f56(g), which can be associated with a 
triangle-type diagram, the time-ordered diagram approach used here gives a factor g\ 
instead of qa as directly obtained from the Weinberg- Tomozawa term. This discrepancy 
points to the fact that the approach misses some contribution. As already mentioned, 
this will be discussed in more details in Sec. 14.21 when making a comparison with the 
expressions obtained from the covariant formalism. 



4 Relation of the EFT approach to the covariant one 

We here discuss similarities and differences between the expressions obtained from the 
full covariant approach and the EFT one presented in sections [2] and [3] respectively. 
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4.1 Similarities (large-M limit) 



To make a comparison of the EFT and time-ordered-diagram approaches with the co- 
variant one, the first step is to derive expressions in the large-M hmit for the last case. 
Taking this limit in the simplest-minded way for the H and G functions given in Eq. (ITOl) . 
one gets: 



{t'-2ml) (Ti X ^ ^ 




Inserting the above limits in Eqs. ([9]), one finds: 



. 4g ~ 2Jt'-Aml 
9uit')M^^ = K—^=K — 



( i\ r> ^ ^ _ r> ^ t' — 27X11. 



g^M3 M'^x I t'-Aml 



_ 7]- 7-' _ 4*7,2 q J-/ _ 277,2 

= ^ + I '—P^) . (18) 

In obtaining the above expressions, one has taken into account that the integral in Eq. 
([9]) for g5e{t') has a higher 1/M order (l/M^). The case of gi^it') is more complicated 
as individual contributions of l/M^ and order cancel. In any case, we notice that 

the large-M limit does not commute with the large-t' limit (compare with the results 
given in Eq. ( JTTl) ). Differences from the "covariant" results are therefore expected when 
considering short distances, where the last limit is relevant. The insertion of the above 
limits in the dispersion-relation integrals, Eqs. ([5]), allows one to recover the expressions 
of the time-ordered-diagram approach at the lowest order, 1/M^, for interactions V44 and 
^56: 



vMm^oo = J^z"".:.. dt' ^ 



I6V2 TC^M^ J Aral ^ {t' + q^) ' 



For the interaction V56, it is noticed that the two terms involving integrands proportional 
to — (t'— 4m^) and 3(t'— 4m^) + 4m^ in Eq. f|T5|) (together with Eq. f|T6l) ) combine to give 
the factor proportional to (t' — 2m^) appearing in the large- M-limit expression of g^eit'), 
Eq. m- 

Expressions for both V44 and V^q can be cast into a form that facilitates the comparison 
with Zhu et a/.'s work [6]. In this order, a factor glj^j^/^AngAM)^ , which can be identified 
as the factor in their work, is partly factored out. Infinities present in the dispersion 
integrals are temporarily ascribed to LEC's, knowing that these ones should be finite in 
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practice. We thus have: 



with: 



44 ' 



V5GiQ)M^oo = V^^iq)+C'^ 



56 ; 



-'56 



(g) = -v^^ (4;;;^M)3 (^?i^(^)-^i(3^(^)-^(^))j- (20) 



The only significant discrepancy with Zhu et al.^s work, Eqs. (jT2ll . concerns the first term 
of V56, which contains a factor g\ instead of qa, confirming the observation aheady made 
in the time-ordered-diagram approach. 



4.2 Differences 

After having shown how the expressions of the EFT (or the time-ordered-diagram) ap- 
proach can be obtained from the covariant ones for the gross features, we now examine 
the differences. 

The first difference concerns the convergence properties of the integral expressions for 
the potentials Vij{q), Eqs. ([5]). While the EFT ones do not converge (infinite LEC's), those 
components produced by the crossed-box diagram, V44 and V34 in the covariant approach 
always converge. This also holds for the crossed-box part of the other components V56 
and V75 but, in these cases, one has to consider a further contribution from the non- 
crossed box diagram. The convergence crucially depends on the way the iterated OPE is 
calculated but there is one choice, quite natural actually, which provides convergence as 
good as for the crossed-box diagram. Though it does not really make sense physically to 
integrate dispersion integrals over t' up to cxd, expressions so obtained provide a reliable 
benchmark, as far as the same physics is implied. 

The second difference concerns the number of components. The covariant approach 
involves many more than the EFT one at NNLO (6 instead of 2). The extra ones imply 
some recoil effect and have a non-local character. They are of higher order in a 1/M 
expansion but, instead, in the large- 1' limit, they compare to the other ones, see Eqs. 

The third difference has to do with the large-M limit of the covariant approach, Eqs. (Q, 
which allows one to recover the structure of the EFT results. The way this limit is taken 
in the H or G functions, or in the factor multiplying the first quantity, is quite rough. It 
assumes approximations like x"^ +4M'^q'^ ~ 4M^g^ {2qT^M > x). Actually, due to the small 
but finite value of the pion mass, there is a very little range of t' values (f— 4m^ < m^/M^) 
where this approximation is not valid. The correction, which disappears in the chiral limit 
(zero pion mass), could affect the long-range part of the interaction. 

The fourth difference involves chiral symmetry and related properties. Contrary to 
what is sometimes thought, fulfilling these properties in calculating the TPE contribution 
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in the covariant approach is possible. This however supposes some elaboration, requiring 
that a description of the NN vrvr transition amplitude entering the dispersion relations 
is consistent with chiral symmetry. In the instance of the Paris model for the NN strong 
interaction this amplitude could be related to experimental data. For the A^A^ weak 
interaction, the strong amplitude was instead modeled from the contribution of a few nu- 
cleon resonances in the s-channel |T3]. This can be essentially achieved by adding to the 
nucleon intermediate state retained here the contribution of the A(1232MeV) resonance. 
This one, in the dispersion-relation formalism, suppresses the low-energy Ntc Ntt 
strong-transition amplitude, otherwise dominated by a well-known large Z-type contri- 
bution inconsistent with chiral symmetry. This part, which involves two pions in an 
isosinglet state (with the cr-meson quantum numbers), is irrelevant here however. Its con- 
tribution is suppressed, in accordance with the Barton theorem [26j which states that the 
exchange of scalar and pseudo-scalar neutral mesons does not contribute to the PV NN 
interaction (assuming CP conservation). This feature largely explains why the PV TPEP 
contribution has not been found as important as originally expected, on the basis of the 
strong-interaction case [TB] . There is another part that is of interest here. It decreases the 
Z-type contribution (the first term of V^qIq) in Eq. (|T5|l ) by an amount which corresponds 
to changing the factor g\ into qa in the first term of V56{q) in Eq. fl20l) . This can be 
checked in the simplest non-relativistic quark model according to the relation: 



,2 



which shows that the contribution of nucleon resonances to the vrA^ scattering amplitude 
can be accounted for by dividing the intermediate nucleon contribution by a factor g\. A 
somewhat different but better argument is based on the Adler-Weisberger sum rule [27] . 
This one, which does not involve any non-relativistic limit, can be cast into the form 1 — 
/ ■ ■ ■ = ^/ g\i where the integral involves the off-mass-shell pion-proton total cross sections. 
For simplicity, we did not retain here the above A contribution, possibly improved for 
other resonances. We nevertheless keep in mind from the previous considerations that the 
discrepancy between the EFT and the covariant approaches, which was noticed for the 
contribution of the triangle diagram in the former one (a factor g^ instead of (7^), could 
be removed by completing the latter one. The corresponding contribution, considered in 
Ref. |T5], amounts to (10—20)% of the one retained here. 

A last remark concerns the comparison of the TPE with the p-meson exchange. For 
a part, the first one was discarded in the past due to possible double counting with the 
second one. With this respect, we notice that the ratio of the V44 and V34 components, 
which could contribute to PV effects in pp scattering {}Sq — ^Pq transition amphtude), 
is very much like the ratio of the local and non-local parts of a standard p-exchange 
contribution. There is some relationship between this result and the fact that the pion 
cloud produces a contribution to the anomalous magnetic moment of the nucleon, which 
compares to the physical one. The problem is different for the other PV transition, ^5*1 — 
■^Pi, where the V56 component could contribute. The corresponding charged-p exchange 
of interest in this case is governed by the PV coupling h'^^ which was predicted to vanish 
in the DDH work [2B]. A small value (-0.7 x lO"'^) was obtained by Holstein ^ by 
considering a pole model. A larger value (— (2— 3)xl0^'') was obtained later on by Kaiser 
and Meissner [30] using a soliton model. In any case, these values lead to negligible 
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effects. We however observe that the PV irNN couphng constant hi. is also small in the 
same models, most of its possible larger predicted value being due to the contribution 
of strange quarks [31]. Sizable values of h'^ could thus be expected. On the basis of a 
dynamical model considering the p meson as a two-pion resonance, one cannot exclude 
values of h'^ in the range of {5—10)^/2 h^., a relation that the results of Kaiser and Meissner 
roughly verify. Contrary to the previous values, the last ones could lead to some double 
counting and it is likely that the h'^ contribution to PV effects is then largely accounted 
for by the two-pion exchange considered in this work. It is conceivable that a similar 
conclusion holds for the other isovector coupling, h^, which contributes to the PV effects 
in pp scattering discussed above. Another aspect of the comparison with a p exchange 
concerns the range of the TPE, which was presented in Ref. ^ as a medium one. This 
could apply to the exchange of two pions in a S" wave {a meson), which contributes to the 
strong A^A^ interaction but is absent in the weak already mentioned. The TPE of 

interest here involves two pions in a P wave with the quantum numbers of the p meson. 
Due to a centrifugal barrier factor, the TPE contribution is shifted to values of t' higher 
than for a S wave, making the range of its contribution closer to a p-exchange one. Some 
numerical illustration is given in next section. 



5 Numerical comparison of potentials 




Figure 4: Spectral functions g{t') in units of GeV~^, with the coefficient K factored out, 
and represented as a function of t'^/^ to better emphasize the low-t' range. Left panel: 
for t'^/^ from threshold to 5 GeV for all functions, to show the relative importance of the 
various components and the beginning of the onset of the asymptotic behavior, 
Right panel: from threshold to 1 GeV for functions gM{t') and gz(,{t') together with their 
large-M limits, to check the validity of this approximation. 

We consider in this section various numerical aspects of the potentials presented in 
the previous one. They successively concern the spectral functions, gij{t'), the poten- 
tials in momentum space, Vij{q), and the potentials in configuration space, Vij{r). In 
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most cases, we directly compare the results of the covariant approach (GOV) with those 
obtained from it in the large-M limit (LM), Eqs. (!20l) . This comparison is more mean- 
ingful than the one with the EFT potential (EFT), Eqs. (fT2!) . as it is not biased by 
the choice of the factor and by the difference of a factor g^, instead of g\, in part 
of the contribution to V^q, of which origin is understood in any case. Before enter- 
ing into details, we notice that the various components of the interaction have a lo- 
cal character for some of them (V44 and V^q) and a non-local one for the other ones 
(V34 and V75). At low energy, however, it turns out that one of the contributions in- 
volving the factor p or p' in their expression, Eq. ([7]), is small. Moreover, with our 
conventions, the spin-isospin factors give the same values for the z component. The 
various components can then be usefully compared, independently of their local or non- 
local character, that is what we do here. Numerical results assume the following values: 
g^j^j^ = 13.45, qa = 1.2695, M = 938.919MeV, = 138.039MeV, rup = 771.1MeV. 




q (GeV) q (GeV) 



Figure 5: Potentials ^44(5), ^'34(5'), v^Q^q) and ^75(5) (together with the OPE one divided 
by 10) for q ranging from to 1 GeV (left panel) and their EFT counterpart for f 44(g) 
and v^elq) (right panel). Ingredients entering the EFT results are specified in the text. 
Notice that the EFT and "covariant" versions of a given component of the potential have 
opposite signs. 

We begin with the first four spectral functions Qijit') entering potentials V44, V34, V56 
and ^75. Their t' dependence is shown for a range of t'^^"^ going from threshold to about 
5 GeV in Fig. H] (left panel). It is aimed to roughly evidence the relative weight of various 
components at small and high t' . At first sight, different potentials have comparable 
sizes. The small-t' range is physically more relevant in the sense that the regime beyond 
1 GeV^ is expected to involve the contribution of other multi-meson exchanges. The 
higher-t' range is more appropriate to illustrate convergence properties. At low values 
of t', one can see some significant differences as expected from the 1/M expansion, Eq. 
( IT8|) . For t'^/"^ < 0.5 GeV, the spectral functions entering the local potentials, V44 and V56, 
dominate those of the non-local ones, V34 and V75. All of them increase in the lower-t' 
range (except in a very small t' range for Vsg) and one has to go to much higher values of 
this variable to observe some saturation and ultimately some decrease. The maximum is 
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Figure 6: Subtracted potentials ^44(5) and ^56(0') for q ranging from to 1 GeV: com- 
parison of the large-M limit approach (LM, Eqs. (1201) ) with the covariant one. The 
asymptotic q dependence of the "covariant" results is a constant one while the one for 
large-M limit results has an extra ln(g) dependence. 



roughly reached around t'^^^ = 2M, which, apart from the pion mass, is the only quantity 
entering the calculations. The decrease, roughly given by up to log terms, ensures 

good convergence properties for potentials Vij{q), Eq. ([5]) (this would not be the case 
for the other Green function mentioned in the text, see details in Appendix IA.2p . The 
validity of the 1/M expansion for the spectral functions 5'44(t') and g^%{t'), which allows 
one to recover the EFT results for the essential part, can be checked by examining Fig. 
m (right panel), where the "covariant" and approximate results are shown for a t' range 
extending to 1 GeV^. It is observed that the dominant term in the 1/M expansion tends 
to overestimate the more complete results both at very low and high values of t' . In 
the first case, the threshold behavior (g^ for g^iit'), axctgim^ / qj^M) for g^^it')) is missed 
(see observation on the 1/M expansion in the previous section). In the second case, 
the overestimation, which is roughly given by a factor t'/M^, tends to increase with t', 
preventing one from getting convergent results. 

In Fig. [5] (left panel), we show the potentials Vij{q) up to g = 1 GeV (together with the 
OPE one that has a strong dependence on q and is divided by 10 to fit the figure). As 
expected from examination of the spectral functions, the TPE potentials have roughly the 
same size and there is no strong evidence that some of them should be more important 
than other ones. It is also noticed that their decrease in the range q =(0—1) GeV is 
slower than the standard p-exchange one, given by l/(m^ + g^), indicating they roughly 
correspond to a shorter-range interaction. The comparison of the potentials, f 44(13') and 
f56(g), with the corresponding EFT ones, Eqs. ( fT2l) . can be made by looking at Fig. [5] 
(right panel). The choice of A^ in the EFT results (A^ = An gAM / gT,^]^) is suggested by 
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r (fm) r (fm) 



Figure 7: Potentials vu{r), vu{r), v^Q^r) and vj^^r) at small and intermediate distances 
(left and right panels respectively). The results shown in the figure represent the above 
potentials multiplied by a phase factor, (in units of GeV~^), to better emphasize the 
most relevant range for applications. 

the large-M limit of the "covariant" expression, see Eqs. (120|) . to make the comparison as 
meaningful as possible. A striking feature appears here: the EFT potentials have a sign 
opposite to the "covariant" ones but the q dependence is roughly the same. This feature 
suggests that the LEG part could play an important role. 

To make a more significant comparison, we subtracted a constant from both potentials 
so that they vanish at = 0. This procedure amounts to using subtracted dispersion 
relations, which leads to convergence in all cases. Moreover, in configuration space, this 
part of the interaction determines the long-range component of the potential, which is 
physically the most relevant one. Results are shown in Fig. [61 The "covariant" and LM 
results have now the same sign. It is however noticed that the present LM results tend to 
overestimate the "covariant" ones. In the limit g — 0, the overestimate reaches a factor 
1.6 for V44 and a factor 1.3 for Vse- This points in this case to the role of higher 1/M- 
order corrections. Actually, the discrepancy vanishes in the limit rriT^/M 0, showing 
that the non-zero pion mass has some effect. The discrepancy tends to slowly increase 
with q (factors 1.9 and 1.5 at q =1 GeV for V44 and V56 respectively), pointing out to 
the role of ln(g) corrections appearing in the large-M limit. Altogether, these results 
are in accordance with the overestimate already noticed for spectral functions in this 
approximation. /^From a different viewpoint, these results confirm the expectation that 
the discrepancy between the EFT and "covariant" approaches shown in the right panel of 
Fig. [5] can be ascribed to contact terms. A rough agreement would be obtained with the 
effective potential obtained in the MS together with a dimensional-regularization scale 
fi ranging from Sm^r to 6171.,^, depending on how this is made (see the definition of this 
scheme at the end of Appendix iBl). 

The Fourier transforms of the Vij{q) quantities, Vij{r), are shown in Fig. [7] for small dis- 
tances (left panel) as well as intermediate distances (right panel). We call them potentials 
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Figure 8: Potentials W44(r) and V5Q{r) at small and intermediate distances (left and right 
panels respectively): comparison of the "covariant" calculations with their large-M limits; 
other definitions or comments as in Fig. [71 The curves corresponding to the large-M limit 
in the left panel tend to 00 when r ^ and have a zero-range part, not shown in the 
figure, with an opposite sign. 

though they are dimensionless quantities, the energy dimension being given by the extra 
operators p. They are multiplied by a factor to emphasize the range which is relevant 
in practice for calculations. At small distances, the comparison of various components 
roughly reflects the one for spectral functions or potentials in momentum space. Exami- 
nation of these results at large distances evidences some significant differences. They have 
a better agreement with expectations from the 1/M expansion or from the very-low-t' be- 
havior of spectral functions. Thus, the local potentials, f44(r) and v^Q^r), have a range 
larger than the other two components, W34(r) and f75(r), do. The differences appear only 
in the range where potentials have small contributions to PV effects. 

The comparison with the LM potentials is given in Fig. [8] for the local components (V44 
and Vse). It is noticed that these last potentials should be completed by contact terms, 
which have a sign opposite to the corresponding curves. Due to the difficulty of represent- 
ing these terms in a simple way, they have not been drawn in this figure. Moreover, they 
are not distinguishable from the LEC's contributions and have thus an arbitrary character 
(they depend on the subtraction scheme). Considering first potentials at intermediate (or 
long) distances, where they can be the most reliably determined, it is found that the LM 
and "covariant" potentials have the same sign despite they have opposite sign in momen- 
tum space. This result is in complete accordance with the subtracted potentials shown 
in Fig. [HI Indeed, the slope of the corresponding results is, up to a minus sign, a direct 
measure of the square radius of the potential weighted by its strength. The negative 
slope for the potential V44 thus indicates that its configuration-space representation is 
positive at intermediate distances (the opposite for V^q). Quantitatively, the significant 
dominance of the LM results over the "covariant" ones (factors 1.7 and 1.4 for V44 and V56 
respectively at r = 0.8 fm) confirms what is found for subtracted potentials in momentum 
space. Considering now the very short-range domain, it is found that the product of the 
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LM potentials with r^, which are shown in Fig. [8] (left panel), diverge like 1/r when 
r ^ 0, up to log factors. The contribution of this part alone to the plane- wave Born am- 
plitude is thus logarithmically divergent. This contribution turns out to be canceled by 
the zero-range one, mentioned above, so that the sum is finite. Thus, there is no principle 
difficulty with this peculiar behavior of the LM potentials in configuration space but, of 
course, some care is required in estimating their contribution. 




r (fm) r (fm) 



Figure 9: Potentials vu{r) and Vr,Q(r) at small and intermediate distances (left and right 
panels respectively): comparison of the "covariant" calculations with a standard rho- 
exchange potential normalized to the same volume integral; other definitions or comments 
as in Fig. [3 

Naively, it could be thought that the TPE is a medium-range interaction, as already 
mentioned. In Fig. [9l we compare potentials f44(r) and f56('") to a standard rho-exchange 
one normalized so that they have the same volume integral. This quantity determines the 
low-energy plane- wave Born amplitude (up to a factor p which can be factored out). At 
very large distances, the effect of the longer-range TPE tail is evident but this occurs in a 
domain where the potential is quite small and will not contribute much. At intermediate or 
even at small distances however, the TPE roughly compares to the p exchange. Actually, 
it turns out to have a shorter range. This reflects the fact that the two-pion continuum 
has an unlimited mass (the integration over t' extends to inflnity). Moreover, it is slightly 
more singular at very small distances, as a result of the extra ln(r) dependence of the TPE 
potential. This last effect is typical of relativistic effects. Both long- and short- distance 
effects can be traced back to the g{t') function which, in the TPE case, extends to both 
small and large values of t' with a maximum around t' = 4M^, while it is concentrated 
around t' = m^ for the p-exchange one ((5(t'— m^) function in the zero-width limit). 

For simplicity, we did not consider explicitly the contribution due to nucleon resonances 
in the two-pion box diagrams, which was accounted for in the 70's in a "covariant" ap- 
proach |T5] or recently in a EFT [T6l [T7] or a TO one [18]. As already mentioned, part of 
it is contained in the EFT approach by relying on the Weinberg- Tomozawa description of 
the TiN scattering amplitude. It only contributes to the V56 component and could repre- 
sent 10—20% of the total contribution, depending on the range and how it is estimated. 
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On top of it, there are further contributions which affect both V44 and V56 components. In 
the case of V44, they amount to an extra 30—40% contribution around 1 fm in the "covari- 
ant" calculation [15] and roughly twice as much in the EFT approach fl6[ [T7] or the TO 
one [12]. For the case of V^q, it is more complicated as the above mentioned contribution 
relative to the description of the vrA^ scattering amplitude has to be disentangled first for 
the "covariant" calculation. When this is done, the extra contribution due to resonances 
decreases and could represent 30—40% of the contribution with nucleons only [15]. This 
is slightly less than what is obtained in the EFT approach [16] . Taking into account that 
the EFT results overestimate the "covariant" ones for the nucleon intermediate state, it 
thus appears that the overestimate for the resonance intermediate state is significantly 
larger. This feature points to corrections of order p^/(M(Ma — M)) affecting the reso- 
nance propagator (dispersion effects), which are known to be important, while corrections 
in relation with the large-M limit are of order p^/M^. It reinforces the conclusion of Ref. 
[T6] that, contrary to the strong-interaction case, the role of resonances, especially the 
A(1232 MeV) one, plays a negligible role in the PV A^A^ interaction. 



6 Estimates of PV effects in two processes 

The TPE potentials considered here have an isovector character. At low energy, they can 
contribute to two different transitions, ^Sq — ^Po, which involves identical particles like 
two protons or two neutrons and ^Si—^Pi, which involves different particles, proton and 
neutron. The interactions V44 and V34 contribute in the first case while the interactions 
V56 and V75 contribute in the other one. Two processes of current interest, where the 
TPE potentials matter, are respectively proton-proton scattering and radiative neutron- 
proton capture at thermal energy. In the first helicity dependence of the cross 
section, Al{E), has been measured at different energies [321 E21 El]- In the second case, 
an asymmetry in the direction of the photon emission with respect to the neutron po- 
larization, A^, has been looked for at LANSCE [35] (the experiment is now running at 
SNS). The TPE contribution to these effects is discussed in the following. The calcula- 
tions have been performed with the A^A^-strong- interaction model, Avl8 [36j, which is 
local wave by wave. Vertex form factors are ignored. On the one hand, the dispersion- 
relation formalism assumes on-mass-shell particles, the contribution due to form factors 
in other approaches being generated by what is included in the dispersion relations. On 
the other hand, the role of form factors was already examined within the EFT approach 
[Tin E]) partly with the motivation of regularizing a potential that is badly behaved at 
short distances, oc {r^^ — c6{f)) where c is infinite and "determined" so that the integral 
over r has a well-defined value. There is no principle difficulty to work with this potential, 
however]^ and we will therefore use it here. This will facilitate the comparison with the 
"covariant" results. 

As a side remark, we notice that form factors different from those mentioned above have 
been used in the context of applying effective field theories to the strong A^A^ interaction 
[1]. They involve a separable dependence of the relative momentum in the initial and final 

^The trick is to separate in the integrands a part determined by wave functions at the origin, of which 
the integral over f is known, the remaining part being well behaved at the origin. 
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states, p and p instead of g = p—p'. Their effect is to smooth out wave functions at short 
distances in accordance with the idea that the corresponding physics, partly unknown, 
should be integrated out and accounted for by LEC's. For such form factors, it would be 
more convenient to work in momentum space. However, in the case where configuration 
space is chosen, the methods we used for dealing with the badly-behaved potential could 
be useful there too. 

6.1 Proton-proton scattering 

The first calculation of TPE effects was done by Simonius [3Z] , with the aim to get some 
estimate for a measurement of PV effects in pp scattering (the Cabibbo model then used 
was not contributing to the pp force in its simplest form). Our results, obtained here 
for three energies at which the longitudinal asymmetry Al[E) has been measured, are 
presented in Table [H 

Examining the "covariant" results, it is found that the contribution of the local term, 
V44(COV), dominates over the non-local one, V34(COV), which appears at the next order 
in the 1/M expansion. The result could be guessed from looking at Fig. [71 Their ratio is 
of the order of the factor 1 + = 4.706 which appears in the p-exchange potential. There 
are reasons to think this result is not accidental (see end of Sect. H]). The present results 
compare to the earlier ones [37] as well as with the value Ai(13.6MeV) = — O.l/i^ that 
has been used in analyses of PV effects [31] . The closeness can be attributed to the fact 
that the Avl8 model employed here and the Reid or Hamada- Johnston models previously 
used are local ones and, moreover, evidence a strong short-range repulsion. Significant 
departures could occur, instead, with models like CD-Bonn [38j or some Nijmegen ones 
[39] . which have a non-local character (see Ref. [ID] for a discussion about the role of 
non- locality). 

Table 1: PV asymmetries in pp scattering at three energies, 13.6, 45 and 221 MeV: 
successively for the potentials V44(COV), V34(COV), V44(rho-like) and V44(LM), in units 
of the hi coupling constant. 



Energy (MeV) 


13.6 


45 


221 


V^44(COV) 


-0.092 


-0.154 


0.072 


l^34(COV) 


-0.022 


-0.042 


-0.029 


V44(rho-like) 


-0.110 


-0.196 


0.096 


1^44 (LM) 


-0.150 


-0.252 


0.127 



Due to its long-range component, one could infer that the TPE contribution should be 
enhanced with respect to the p-exchange one when the effect of short-range repulsion is 
accounted for. The comparison of V44(rho-like) and V44(COV) results shows that this is 
the other way round. This can be explained by the fact that the long-range contribution 
where the TPE dominates over the p-exchange has little contribution to the asymmetry 
(a few %). Instead, the short-range contribution where the TPE dominates over the p- 
exchange plays a bigger role. The effect of short-range repulsion on this contribution will 
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therefore be enhanced, hence decreasing of the total TPE contribution with respect to 
the p-exchange one. 

6.2 Asymmetry in neutron-proton radiative capture 

Contrary to pp scattering discussed above, the asymmetry involves a non-zero contri- 
bution from OPE. Thus, independently of the value of the coupling /i^, one can directly 
compare the TPE and OPE contributions. The last one has been extensively studied, see 
Ref. pi] and earlier references therein, and Refs. [TUl IHl |12] for more recent works. Its 
contribution is approximately given by A^f{OPE) = — O.ll/i^ (= — 0.112/i^ for the Avl8 
model used here). 

The TPE contributions to for the "covariant" case are given by: 

A^iV^eiCOV)) = 0.0093/1^, 

A^iVr^iCOV)) = -0.0040/1^, (22) 

while those for the rho-like and large-M limit are: 

A^(l4,6(rho-like)) = 0.0093 /i^, (23) 
A^{V,e{LM)) = 0M41hl. (24) 

Considering the "covariant" results, it is first noticed that the contribution of the local 
term, V56(COV), dominates over the non-local one, V75(COV), which appears at the next 
1/M order. The effect is however less important than in pp scattering, which can be 
inferred from looking at the corresponding potentials in Fig. [7l Moreover, contrary 
to this process, their contributions have opposite signs. Thus, the total contribution 
represents only —5% of the OPE one. This is good news in the sense that it confirms that 
the asymmetry A^ represents the best observable to determine the coupling h]^. Present 
results reasonably compare to an earlier one [19]: A^{TPE) = 0.008 /i^, which is added 
up by A^iV^e) = 0.0107 /i^ and ^^(1^75) = -0.0027 /i^. Part of the difference for A^iV^e) 
can be traced back to the omission here of the A resonance contribution in intermediate 
states. It is reminded that this contribution tends to make the vrA^-scattering amplitude 
consistent with the Weinberg- Tomozawa coupling, resulting in an enhancement of the 
V56 interaction (see discussion in Sec. H]). The main difference for ^4^(1^75) is due to the 
extension of the integration over t' in the dispersion relation from 50 to 00 here. 

At first sight, the comparison of the TPE and p-exchange results, which are essentially 
the same, shows features different from those observed in pp scattering. Examination of 
Fig. [9] indicates that, in comparison to the V44 potential, the enhancement of the TPE 
potential, V56(COV), over the p-exchange one, V56 (rho-like), is larger at large distances 
and smaller at short distances. As a result, the two effects due to the short-range repulsion 
mentioned for pp scattering tend to cancel here. 
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6.3 Results in the large-M limit 



Involving the subtraction of an infinite term, results obtained in the large-M limit, which 
can be identified with the EFT ones for the essential part, cannot be directly compared 
to the "covariant" ones. A first look nevertheless shows that the results, apart from 
an enhancement by a factor 1.5 or so, are very similar. Moreover, from considering the 
plane-wave Born approximation, we should have expected opposite signs. It is therefore 
appropriate to give some explanation helping to understand these results. 

We first checked which range was contributing to PV asymmetries and found that the 
role of the range below 0.4—0.5 fm was rather small (10—20%). The observed enhancement 
of the large-M result over the "covariant" one thus reflects the similar enhancement that 
can be inferred from the corresponding potentials in Fig. [8|, around 0.8 fm. The little 
role of the range below 0.4—0.5 fm is not expected from considering Fig. [3 Instead, it 
can be related to the strong interaction model, Avl8, used to calculate wave functions 
entering estimates of observables. This model evidences the effect of a strong repulsion 
at short distances. As this effect is much less pronounced with non-local models, one 
can expect that the use of models like CD-Bonn or some Nijmegen ones, will show some 
differences with the present results. We can anticipate an increase of the magnitude of 
the "covariant" results and a decrease of the LM ones (without excluding in this case a 
change in sign reminding that one for the plane- wave Born approximation). 

While trying to understand the role of higher 1/M-order corrections, we face the prob- 
lem that the corresponding contributions to potentials are more singular than those for 
V44(LM) or V56(LM), requiring some regularization and introduction of further LEC's. 
Restricting our study to the range r > 0.4—0.5 fm, which provides most of the contribu- 
tion to the PV asymmetries, we found that the contribution at the next 1/M order had 
a sign opposite to the dominant one, which it largely cancels in the range around 0.8 fm. 
The correction, which is larger than needed, suggests that other corrections are neces- 
sary to ensure reasonable convergence. Among them, one could involve chiral symmetry 
breaking. We checked that in the limit of a massless pion, the "covariant" potential f44(r) 
and its large-M limit are significantly closer to each other. A non-zero pion mass could 
thus explain a sizable part of the discrepancy for observables between the "covariant" and 
the large- M-limit results. We notice that the comparison of the EFT and "covariant" ap- 
proaches in the strong interaction case evidences features similar to the above ones (see 
Ref. [21] and references therein). As far as we can see, the similarity is founded for a 
part. 

The large-M limit of "covariant" potentials assumes a particular subtraction scheme. 
In another scheme, the MS one, part of the functions L{q) appearing in Eqs. (l20l) are 
replaced by L(g)— 1— ln(/i/m,r) (see Appendix[B]). Due to the effect of short-range repulsion 
in the Avl8 model, the correction has little effect on the calculated observables (a few % 
increase for the S—P transition part). This would be different for non-local models where 
the correction could be significantly larger. Depending on the value of the parameter /x, 
a large part of the sensitivity of the LM results to strong-interaction models mentioned 
above could be removed [10]. We notice that a discussion similar to the above one on 
the role of the subtraction scheme has also been held in the strong- interaction case [1|. 
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It was concluded that the spectral function regularization scheme (SFR) was providing 
better convergence properties than for the dimensional-regularization one, due to avoiding 
spurious short-range contributions. To some extent, the SFR scheme is close to the MS 
one considered here, with the parameter A introduced in the former case being replaced 
by the quantity 2/i in the latter case. 



7 Conclusion 

In the present work, we have compared different approaches for incorporating the two- 
pion-exchange contribution to the PV A^A^ interaction. They include a covariant one, 
which fully converges and can thus be considered as a benchmark, and an effective-field- 
theory and a time-ordered one which can contain infinities. The last two approaches 
involve two components at the leading order, with a local character, while the covariant 
one involves both local and non-local. For a given transition, one can thus compare the 
local components obtained from different approaches on the one hand, local and non-local 
components on the other hand. These two comparisons can allow one to assess how good 
is the assumption of dominant order in the effective-field theory approaches as well as the 
role ascribed to LEC's. 

We first notice that the effective-field theory (EFT) approach, the time-ordered-diagram 
approach and the limit of the covariant one at the lowest non-zero order in the 1/M ex- 
pansion (LM) essentially agree with each other for the local terms. Possible discrepancies 
involve ingredients that have been omitted (contribution of baryon resonances in particu- 
lar) but are unimportant for the comparison. We can thus concentrate on a comparison 
of the covariant approach with its LM limit. Taking into account that this approach is 
determined up to contact terms, rough agreement is found. This is better seen by consid- 
ering the subtracted potentials in momentum space or intermediate distances (r = Ifm) 
in configuration space. Quantitatively, the LM (EFT) approach tends to overestimate 
the "covariant" results. At low q or at intermediate distances, the effect reaches factors 
1.7 and 1.3 respectively for the transitions ^5*0 — ^Pq and ^Si — ^Pi. At very small but 
finite distances, the LM (EFT) potentials become very singular and their contribution to 
physical processes diverge. This divergence is canceled by the contribution associated to 
the contact term so that the total result is finite (after renormalization). This peculiar 
behavior at r = and around contrasts with the smooth but diverging behavior in mo- 
mentum space of the published EFT two-pion-exchange interaction. In this case, it turns 
out that the sign of the potential is opposite to that one at finite distances, which is the 
most relevant part. This suggests that this EFT two-pion-exchange potential is domi- 
nated by an unknown contact term, as far as a comparison with the "covariant" result is 
concerned. The problem disappears with a different subtraction scheme, like the minimal 
one with a dimensional-regularization scale /i in the range (3—6) m^. This last choice 
tends to minimize the role of short distances, confirming the absence of a large sensitivity 
to cutoffs observed elsewhere [10] . Interestingly, the MS scheme corresponds to cutting off 
the dispersion integrals at a value of t'^/^ around 2/i, which, together with the above value 
of /X, is about (1—2) GeV. This is quite a reasonable value for separating the contribution 
of known physics from the unknown one to be integrated out. We thus believe that the 
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choice of the MS scheme would be more appropriate, the interaction then ascribed to the 
EFT two-pion-exchange one being a better representation of the most rehable part of the 
two-pion-exchange physics, which occurs at intermediate and large distances. 

Comparing the non-local components to the local ones, it is found that they are rather 
suppressed at large distances. For such distances, the contribution to dispersion relations 
is expected to come from values of t' smaller than M^. The suppression of non-local terms 
then reflects the fact that they have an extra 1/M factor in the 1/M expansion. At short 
distances, instead, the contribution to dispersion relations comes from large values of t' . 
In this case, the local and non-local terms have the same 1/M order and they tend to 
have comparable contributions. Another aspect of the dependence on the range concerns 
the comparison with the p-meson exchange. Not surprisingly, the two-pion exchange 
dominates the p-meson exchange one at long distances but, the potential being relatively 
small there, not much effect is expected from this part on the calculation of observables. 
The two-pion-exchange contribution is slightly dominated by the p-meson exchange one 
at medium distances and dominates again the last one at very short distances. Taking 
into account that the dominant contributions come from short and intermediate distances, 
the two-pion-exchange contribution turns out to have a shorter range than the p-meson 
exchange. 

We looked at the TPE contribution in two physical processes, pp scattering and radiative 
thermal neutron-proton capture. Roughly, they confirm what could be inferred from 
examining potentials. The results from the "covariant" approach, calculated with the 
Avl8 A^A^ strong-interaction model, essentially agree with earlier estimates based on other 
models. The main discrepancies evidence the role of some inputs such as the restriction 
on the t' value in the dispersion relations or the role of resonances in modeling the ttA" 
scattering amplitude which enters these relations and we omitted here for simplicity. 
Despite plane-wave Born amplitudes calculated with the EFT two-pion-exchange and 
the "covariant" approaches have opposite signs, it turns out that their contributions to 
observables are relatively close to each other. This feature points to the Avl8 model, 
which produces wave functions that evidence the effect of a strong repulsion at short 
distances. Such a property is interesting in that the main contribution to observables 
comes from intermediate and long distances, where the derivation of the TPE potential 
is the most reliable. It is thus found that the EFT two-pion-exchange at NNLO tends 
to overestimate the "covariant" results by about 50%, pointing to a non-negligible role 
of next order corrections. This is confirmed by the consideration of non-local terms, 
which correspond to higher-order terms. Their contribution is especially important in 
neutron-proton capture, due to a destructive interference with the dominant one. Thus, 
the result at the dominant 1/M order (LM) in this process exceeds the "covariant" one 
by a factor of about 2—3. Interestingly, present results are rather insensitive to the 
subtraction scheme as far as the coefficient, l-|-ln(p/m,r), remains in the range of a few 
units. This is a consequence of the strong short-range repulsion present in the Avl8 
model. We however stress the fortunate character of this result. The use of non-local 
strong-interaction models, like CD-Bonn or some Nijmegen ones, could lead to different 
conclusions. Actually, far to be a problem, the dependences on the model expected for the 
contribution of the term l-|-ln(p/m7r) and the EFT potential calculated in the maximal- 
subtraction (MX) scheme are likely to cancel for a large part. 
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By comparing different approaches to the description of the PV TPE A^A^ interaction, 
it was expected one could learn about their respective relevance. Implying a natural cut- 
off of the order of the nucleon mass, the "covariant" approach provides an unavoidable 
benchmark. This is of interest for the EFT approach which, up to now, has been con- 
sidered at the lowest 1/M order. In improving this approach, a first step concerns the 
subtraction scheme. The minimal-subtraction one (MS), which involves a renormaliza- 
tion scale, /i, is probably more favorable. By taking for this scale a value of the order of 
the nucleon mass (or the chiral-symmetry-breaking scale A^^), the scheme better matches 
the separation of the interaction into known and less known contributions, corresponding 
respectively to long and short distances. The LEC's so obtained could be less dependent 
on the strong-interaction model. The next step should concern higher 1/M-order terms, 
whose contributions are not negligible. This is likely to require a lot of work and cau- 
tion, as the singular behavior of these terms at short distances increases with their order. 
Meanwhile, the "covariant" results could provide both a useful estimate and a relevant 
guide for their study. 
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A Subtraction of the iterated OPE and related ques- 
tions 

A.l Expressions of the spectral functions, g{t') 

Historically, the derivation of the isovector PV two-pion exchange started with the cal- 
culation of the crossed diagram [12]. The calculation of the non-crossed diagram, which 
implies the removing of the iterated OPE contribution and thus requires more care, came 
slightly later [ISjEl Last works along the same lines [HI [15] considered the two types of 
diagrams on the same footing. We first remind here some results relative to the crossed 
and non-crossed diagrams with the notations of Ref. [13] (functions gAit'), 9 sit'), 9c(t') 
and goit'), QEit') respectively). The functions QAit') and goit') correspond to the same 
spin-isospin structure. The two versions of the iterated OPE discussed in the text, which 
concern the Qoit') and gE{t') functions, are considered. The one employed in Ref. [T3] 

^The paper contains editor mistakes that could obscure its understanding: the contents of figures 1 
and 2 should be interchanged and the number -1.61 in the table should be replaced by -0.61. 
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corresponds to the term with the factor {E + M) in the integrand while the other one 
considered in later works contains the factor 2E. The expressions read: 



9Ait') 



1 fG Hx ^ 



9 Bit') = -^9A{t'), 




, , . X f nii \ G 



M2m2 ^V2Mg^/ 2xM 
k\ dk^ f 2E oi (E+M) ^k'^-x E-M^^ x 



p,/pt'-(m2+A;2)2 V ^ 4M 2 ^ 2M2 



r^l dk^ 



+ 



k^JkH'-iml + k'^) 



2 



( 2E oi {E+M) Ak'^-xf 2{E-M){k^-x) ^{E-Mf^ 2x' 



\ E^ V Mt' t' M > MH 



(25) 



where the various functions, q.,,, x, Xi G, H, are given in the text, Eq. ( ITOl) . The 
writing slightly differs from Ref. [IBj. No non-relativistic approximation is made for 
the integrands. Moreover, the original term, arctg(2Mg^/m2)/m2, has been transformed 
into (71/2 — arctg(m2/2Mg7r))/'^^ and the factor Ti/2/m1 has been inserted in the integral 



using the relation J^^^ dk'^/yk'^yk'^t' — {m'^ + k'^y j = Ti'/m'^. With this rearrangement, it 
can be checked that the integrand has no singularity at /c^ = 0. 

The g{t') functions considered in the present work are related to the above ones by the 
relations: 

g^^{t') = Kgcit'), 
gsiit') = Kg Bit'), 
g,,{t') = K (gAif) + gnit')) , 

gj,{t') = KgEit'), (26) 



where K is an overall constant given in Eq. ffTOj) . An important point to note is that 
the contributions of the term G/2xM in gA{t') and gnit'), which dominate at low energy, 
exactly cancel. This cancellation is important in restoring the crossing symmetry for 
pions, a property that is fulfilled by the effective pion-nucleon interaction introduced in 
the EFT approach (triangle diagram). 
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A. 2 Asymptotic behavior of the g{t') functions 



The asymptotic behavior of the spectral functions g4,i{t') and (734 (t') 5 which only involve 
the crossed- diagram contribution, can be easily obtained from their expressions, Eqs. ([9]). 
The dominant term is of the order l/{M\/t/) up to some log factors, see Eq. (fTTI) . The 
asymptotic behavior of the two other spectral functions, g^%{t') and gj^it'), is considerably 
more complicated. Their analytic part contains terms with the behavior —Ky/t' /2M^ and 
Ky/f/M^ respectively while the integral part requires careful examination. 

The dominant term in the integral is given by the part proportional to x and in the 
integrands of g^eit') and gi^it') respectively. As the integrands are the same up to a factor 
— 4x/t', it is sufficient to consider the first case. Its contribution becomes: 



'56 



-Kx 



k'i dk"^ ( 2E {01 E + M) 1 



k^^kH'-{ml + k^Y\ 2MV 

~ r^l dk^ f{E-M)^ {E-M){2E + M)^ 

k^JkH'-iml + k^)A'^M^y''' UPY^ - 



K,f''^ ( I fo, (2^ + ^) )] (27) 

-L2 /p^/_(^2^p)2 \2M^E{E + M)^ iM^ E^{E + M)^ J ' ^ ' 



k1 

where the first case corresponds to the quadratic-energy dependent Green's function re- 
tained here and the second case to the linear one. After some algebra, it is found that in 
the large-t' limit and a negligible pion mass, the integral with the first integrand writes: 

^ 00) . ^ 00) . - + i - il,>(2)) . (28) 

The absence of the intermediate term 1/(M^) results from a non-trivial cancellation. All 
the other terms not considered here behave like l/(M-\/t') (up to log terms) for the most 
important ones. It can thus be checked that the contributions of order \/V/M^ to the 
spectral functions g^Q{t') and g75{t') cancel, leaving contributions of order l/{M^/t'). 

Have we used the other Green's function, the dominant contribution resulting from 
performing the integral would read: 

h,{t' ^oo)c^ -hj,it' - 00) ^ ^ , (29) 

while the corresponding g{t') functions would be given by: 



g,,{t' ^00) :^ -^gr,{t' 00) ^ . (30) 

Thus, contrary to the quadratic-energy dependent Green's function, no cancellation of the 
dominant terms is found with the consequence that the dispersion integrals involving the 
spectral functions g^eit') and gj^it'), Eqs. ([5]), do not converge (logarithmic divergence). 
Moreover, these functions evidence a sign change around t' =(20—25) GeV^. Actually, this 
has little effect on the spectral functions at low values of t' but, of course, this prevents 
one from getting convergent results. The sensitivity of the configuration-space potentials 
was not exceeding 10% at r = 1 fm [T5] . 
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B PV two-pion exchange NN potential from EFT 



In this Appendix, we give the detail of the momentum-space contributions from the two- 
pion exchange diagrams shown in Fig. [21 employing heavy-baryon chiral Lagrangian and 
the dimensional regularization in the d-dimensional space-time, d = 4 — 2e. 

/^From the diagram (b), (c) and (d), we get: 



V, 



X 



-i{Ti + T2) ((Ti X 0-2) -g 



1 u 
- - 7 + ln(47r) + In 
e V m; 

1 

e 



-i[fi X T2) (ai + (?2)-g 



X 



V, 



id) 



-i{Ti X T2} {(Ti + cr2)-q 



2V27rg\hl 

V27rg\hl 
2(47rA)3 

- 7 + ln(47r) + In 
2(47rA)3 




(31) 



+ 2 - 2L{q) 



ml 



+ 



+ ?,L{q)-H{q) 



(32) 



X 



- - 7 + ln(47r) + In 




+ ^L{q)-H{q) 



(33) 



where L{q) and H[q) are defined in Eq. f|T3|) . and the Euler number is given by 7 = 
0.5772 The quantity represents the scale of the dimensional regularization, and 
q = \q \ with q defined by Eq. ([21). We note that we have employed, in the calculation of 
the diagram (d), the same prescription to subtract the two-nucleon-pole contribution as 
in Ref. ^, Eq. (C.3) . 

All these potentials from the loop diagrams have infinities, i.e., the 1/e terms. Along 
with the finite constant terms, they are renormalized by the counter terms (PV NN 
contact terms) C2 + C4 and Cq (see Eq. (5) in Ref. [6]). In the maximal-subtraction 
scheme (MX), the L{q) and H{q) terms, which are the only ones to contribute to the 
finite-range potential, are retained. This procedure fixes the associated contact term. 
In the minimal-subtraction scheme (MS), besides L(g) and H{q) terms, an extra term 
Ini^jj^/m^ + 2 (rather written as 2(l-|-ln(/i/m^)) in the main text), which produces a 
contact term, appears. 
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